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waters of hydration in mechanically strained elastin 
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We report on the direct measurement of the correlation times of the protein backbone carbons 
and proximal waters of hydration in mechanically strained elastin by nuclear magnetic resonance 
methods. The experimental data indicate a decrease in the correlation times of the carbonyl carbons 
as the strain on the biopolymer is increased. These observations are in good agreement with short 
4ns molecular dynamics simulations of (VPGVG)3, a well studied mimetic peptide of elastin. The 
experimental results also indicate a reduction in the correlation time of proximal waters of hydration 
with increasing strain applied to the elastomer. A simple model is suggested that correlates the 
increase in the motion of proximal waters of hydration to the increase in frequency of libration of 
the protein backbone that develops with increasing strain. Together, the reduction in the protein 
entropy accompanied with the increase in entropy of the proximal waters of hydration with increasing 
strain, support the notion that the source of elasticity is driven by an entropic mechanism arising 
from the change in entropy of the protein backbone. 

PACS numbers: 87.64.kj, 87.80.Lg 



Elastin, the principal protein component of the elas- 
tic fiber, is a remarkable biopolymer that gives rise to 
the resilience of many vertebrate tissues [iHSj. The 
structure-function relationship of this system has been 
of interest since Partridge's original work in isolating 
elastin from other tissue constituents and a model in- 
volving cross-links of hydrophobic and hydrophilic do- 
mains 15 . Central to the elasticity of elastin is the de- 
gree of protein solvation as well as the polarity of the 
solvent as made evident by experimental studies of the 
Young's modulus-elastin only exhibits elasticity when 
hydrated A 'liquid drop elastomer' model was sug- 
gested wherein both the surface interaction between hy- 
drophobic and hydrophilic phases and the change in con- 
figurational entropy contribute to elasticity [5] . The no- 
tion that the elasticity of elastin is entropic in origin was 
established early on and a model was proposed that in- 
corporated a network of random chains [7]. However, 
experimental work gave evidence of ordered structures in 
elastin, specifically, the presence of both a-helices and 
/3-turns [1,8 . The short model peptide (VPGVG)„, an 
abundant motif of elastin, has been extensively studied 
as it exhibits the inverse temperature transition that is 
characteristic of the hydrated protein U [3] • In studies 
of the (VPGVG)„ peptide, a model was proposed that 
incorporates a hydrogen bond formed between Vail and 
Val4 residues, resulting in suspended segments between 
/3-turns [TU]. In this model, termed the Librational En- 
tropy Mechanism (LEM), the suspended segments are 
free to undergo large-amplitude low-frequency torsional 
oscillations; the amplitude might be damped during ex- 
tension, giving rise to a decrease in entropy of the seg- 
ment with increasing strain, providing the driving mecha- 
nism for elasticity [10] . In addition, type II /3-turns have 
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also been proposed for the (GXGGX) sequence, a sec- 
ond repeating motif of elastin, wherein hydrogen bonds 
are formed either between Glyl and Gly4 residues or be- 
tween the second and fifth X residues; the interchange 
between the two conformations results in a dynamical /3- 
turn sliding giving rise to elasticity [2]. In both model 
sequences, according to the LEM model, the elastic force 
is borne by the protein backbone by virtue of a decrease 
in entropy upon extension. In all elastin models, wa- 
ters of hydration arc tacitly assumed and yet the contri- 
bution of proximal waters of hydration to the entropic 
mechanism of elasticity has eluded experimentalists. Re- 
cent short time Molecular Dynamics (MD) simulations 
of (VPGVG)i8 have suggested that the orientational en- 
tropy of waters hydrating hydrophobic groups decreases, 
whereas the elastin backbone is more dynamic and has 
a greater entropy upon extension than in a relaxed state 
in disagreement with the LEM model [9]. In this work 
we report on the first direct measurement of the motional 
correlation times of proximal waters of hydration and the 
protein backbone of mechanically strained elastin. 

Nuclear Magnetic Resonance (NMR) experiments were 
performed on purified bovine nuchal ligament elastin 
purchased from Elastin Products Company, LLC. The 
elastin was purified by the neutral extraction method of 
Partridge [1], and was suspended in a mixture of D2O 
and H2O with a volume ratio of 50:50 with 0.0003g/mL 
of sodium azide added as a biocide. The elastin samples 
were mechanically strained in 1.5cm long, 5mm diameter 
NMR tubes with both ends held by polytetrafluoroethy- 
lene tape and sealed using ethylene-vinyl acetate. Four 
samples were prepared; the strain (AL/L) on the samples 
was 0, 17, 35 and 43 percent, and are denoted I, II, HI, 
and IV respectively in the remainder of this work. Addi- 
tionally, after the experiments were performed on sample 
IV, the sample was released and measured after 72 hours 
and is denoted sample V. The loss of water in any of the 
samples was less than 1 percent over the entire course of 
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experiments. All the experiments were carried out on a 
200MHz Tecmag Apollo NMR spectrometer using a Doty 
solids NMR probe. The 90° pulse time was 9/xs in the 
^■^C measurements and 19/iS in the measurements. In 
the ^'^C measurements, a radio frequency field strength 
of We =27. 7kHz was applied with SPINAL-64 decou- 
pling during signal acquisition '11]. All the experiments 
were conducted at approximately 25° C. 

We employed the 2D Ti-Tz NMR technique [12], and 
measured the Ti and T2 relaxation times of water 
in hydrated elastin. By applying a 2D Inverse Laplace 
Transform (ILT) of the experimental data, the Ti and T2 
relaxation times are correlated and are manifested into 
a peak in the 2D map. The 2D ILT map of sample I is 
shown in Fig. 1 as an example. Four peaks are observed 
in the figure, corresponding to four reservoirs of water 
with distinguishable dynamical properties; the 2D ILT 
maps for all samples studied showed four components. 
For a spin 1=1 nucleus, such as ^H, the relaxation is 
governed by the quadrupolar interaction and the Ti and 
T2 are given by 

^ = ^(1 + '^mJic.n) + 4J{2u^d)} (1) 



^ = ^(1 + j)C,{9J{0) + 15J{u;d) + 6J{2u:d)} (2) 

where ujjj is the Larmor frequency for ^H, Cq = 
and 77 is the asymmetry parameter of the potential V and 
assumed to in our case. The spectral density is given 
by J{io) = i+(I)'t )^ with Tc defined as the correlation 
time of the fiuctuating quadrupolar field of the nu- 
cleus 13 . Using Eq. 1 and 2 and our measured values of 
Ti and T2 , the Tc was determined for the four components 
shown. Component ai is free water, by virtue of the fact 
that Ti ~ T2 , a2 is water that resides between the elastin 
fibers and /3 is a mobile component that resides within 
the elastin fiber as verified by a T2-T2 exchange experi- 
ment ^14j . The water corresponding to the peak with the 
shortest Ti and T2 relaxation times, labeled 7 in Fig. 1, 
has been determined to be in closest proximity to elastin, 
also by the T2-T2 exchange experiment and by consider- 
ing the correlation times of the slowest two components 
in the 2D ILT map (rf =48.8ns and r2'=111.7ns). In this 
work, we therefore study the changes in the correlation 
time of component 7 as a function of mechanical strain. 
For samples I, II, III and IV, the correlation times are 
plotted in Fig. 2a as a function of AL/L. It is evident 
from the figure that the correlation time of proximal wa- 
ters of hydration decreases with increasing strain on the 
biopolymer. Not shown in the figure is the correlation 
time of the water from sample V, which was measured to 
be the same, within our experimental uncertainty, as that 
measured in sample I. This observation indicates that the 
dynamics of the water in close proximity to the protein 
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FIG. 1. 2D ILT map of the T1-T2 NMR relaxation times 
of water in sample I, defined in the text. Four components 
are discernable and are labeled ai, 0:2, /3 and 7. The dashed 
lines are used to guide the eye for the region of the 2D map 
where Ti is approximately equal to T2. The signal intensity, 
indicated by the colorbar, is shown on a logarithmic scale. 



I- 































5 


10 


15 20 25 

AL/L [%] 


30 


35 


40 45 


^- 




























I - 




















---m 






i 



5 10 15 20 25 30 35 40 45 
AL/L [%] 



FIG. 2. Correlation times measured in samples I, II, III and 
IV as defined in the text, (a) Correlation time of the fiuctuat- 
ing quadrupolar field of D2O in close proximity to elastin, as 
determined from Eq. 1 and 2 and the measured '^H Ti and 12. 
(b) Correlation times characterizing the fiuctuating dipolar 
field experienced by the carbonyl (open circles) and aliphatic 
carbons (closed circles) of elastin, as determined from Eq. 3 
and 4 and the measured ^^C Ti and Tip. The dashed lines 
are used to guide the eye. 



returned to what was observed in an unstrained sample 
when the applied strain was removed. 

Natural abundance ^^C Direct Polarization (DP) and 
— ^^C Cross Polarization (CP) NMR experiments 
were performed |15| , and exemplary spectra are shown 
in Fig. 3 from sample I. A phase cycling scheme was 
implemented in the CP experiment, such that the result- 
ing signal observed on the ^^C spectra results only from 
cross-polarized magnetization. Two features are evident 
in the spectra. First, it is clear from the DP spectrum, 
shown in Fig. 3a, that the signal for the carbonyl car- 
bons (centered at 173ppm) is clearly distinguishable from 
that of the aliphatic carbons (16ppm-60ppm) even with- 
out any magic angle spinning [TB]; this allows us to dis- 
tinguish these two signals and measure the respective re- 
laxation times unambiguously. Second, the overall signal 
in the CP spectrum (Fig. 3b) is smaller than that of the 
DP spectrum (Fig. 3a) , indicating a high degree of mo- 
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FIG. 3. Natural abundance ^^C NMR spectra of sample 
I defined in the text, (a) Direct Polarization (DP) and (b) 

^^C Cross Polarization (CP). The contact time for 
— >■ '^^C CP was 1.4ms. The data shown were accumulated 
with 14512 scans, with a recycle delay of 5s. All spectra are 
referenced to adamantane. 



bility for this system (a factor of 3.977 in gain is expected 
for cross polarization for the case of a rigid solid [15j). 
Similar spectral resolution was achieved on all samples 
studied. 

Using the DP technique, we obtained the Ti and Tip 
times for the carbonyl and aliphatic carbons of mechani- 
cally strained elastin. Given the low spectral resolution of 
the static sample, one cannot distinguish different signals 
in the aliphatic region and the dynamics of individual 
constituents may vary. However, the ^^C Ti and "'^H Tip 
values measured on hydrated elastin using magic angle 
spinning methods at a similar temperature showed sim- 
ilar values for the different components of the aliphatic 
region [16] . It is therefore reasonable to study the region 
by employing an effective Ti and Tip relaxation time. In 
hydrated elastin, it has been demonstrated that the car- 
bonyl chemical shift anisotropy is averaged to a great 
extent due to the high mobility of the system [16]. The 
^^C relaxation times are thus mediated largely by carbon- 
proton dipolar interactions and the ^■^C Ti and Tip are 
given by 

^ = ^,CD{J{iOH-ioc) + 3J{uJc)+6J{^H+^c)} (3) 
J 1 iU 



Ti 



ip 
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+&J{LOH) + &J{iOH+iOc)} (4) 



where luh and ujc are the Larmor frequencies for and 
respectively, and Co = Clf^ f [11]. In the spec- 
tral density J{uj) in Eq. 3 and 4, Tc is the correlation 
time of the fluctuating dipolar field experienced by the 
^■^C spins, which is reflective of the overall motions of the 
^■^C and nearby nuclei. Using Eq. 3 and 4, and the 
measured relaxation times the Tc is determined for the 
carbonyl carbons and aliphatic carbons in all the sam- 
ples. The results are presented in Fig. 2b as a function 
of the strain on the sample. It is clear from the figure 



that the correlation times for both carbonyl and the over- 
all aliphatic carbons decrease with increasing strain. The 
physics of this result dictates an increase in the fluctua- 
tion of the dipolar field surrounding both carbonyl and 
aliphatic regions with increasing strain. The measured 
correlation times of the carbonyl and overall aliphatic 
carbons of sample V, not shown, were within the experi- 
mental uncertainty of that observed in sample I. To help 
visualize possible dynamical properties of strained elastin 
measured in our experimental work, we have studied the 
elastin mimetic sequence (VPGVG)3 by MD simulation. 

The MD simulations were performed at 25° C using the 
OPLS-AA/L force field model |lSj and a Berendsen ther- 
mostat [IS] in GROMACS [20]. The peptide was termi- 
nated by (-NH-CH3) and (-CO-CH3) groups at its C- and 
N-terminus, respectively. The model system was placed 
in a cubic box solvated with water using the SPC216 wa- 
ter model |21| and was run for 4ns. Two scenarios were 
studied: (a) without an applied pull force and (b) with 
the Glyl4 atoms being fixed and a pull force being ap- 
plied on the atoms of Val4 with a constant force equal to 
50000 kJ/mol nm. The resulting increase in the radius 
of gyration, as determined by the position of the C", is 
highlighted in Table 1 as well as the increase in the num- 
ber of water molecules in 0.6nm of the peptide resulting 
from the applied strain. Further details of the mechan- 
ical properties of this elastin mimetic peptide, and oth- 
ers, will appear in a forthcoming manuscript. For each 
scenario the eigenvalues, Xu, of the mass weighted co- 
variance matrix of the MD trajectory were computed. A 
quasi-harmonic approximation was applied whereby the 
frequencies, uji of uncorrelated simple harmonic oscilla- 
tors are given by u.i — \J^^^ 122] • The entropy of a 

harmonic oscillator with frequency uii at a temperature 
T is given by 



ln{l - e-'^'^''^^)] (5) 



The histogram of the frequencies derived from the simu- 
lations are plotted in Fig. 4 for both scenarios studied. 
The resulting total entropy of the peptide was determined 
by summing the entropy contribution of each frequency 
and is tabulated in Table 1. Fig. 4 shows a shift in pop- 
ulation to higher frequency when the peptide is strained. 
According to Eq. 5 high frequencies contribute less to en- 
tropy, thus Fig. 4 implies a decrease in the total peptide 
entropy when strained; this is observed in our results tab- 
ulated in Table 1. In the LEM model [10], the peptide 
is treated as a series of quasi-harmonic oscillators un- 
dergoing low-frequency, large torsional oscillations when 
relaxed; the amplitude of the motion is reduced upon ex- 
tension. Our experimental data, highlighted in Fig. 2b, 
indicate that there is a decrease in the correlation times 
of both carbonyl and aliphatic carbons, implying that the 
motion of the fluctuating dipolar field of both regions in- 
creases with increasing strain. We attribute the decrease 
in the measured correlation times to an increase in fre- 
quency of the protein backbone motion as observed in our 
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FIG. 4. Histogram of frequencies derived from the quasi- 
harmonic approach in the MD simulations of the elastin 
mimetic peptide (VPGVG)3 under relaxed (blue line) and 
strained (red line) states. 





Relaxed 


Strained 


Peptide entropy [kJ/mol-K] 


1.63 


1.46 


No. of H2O within 0.6nm of peptide 


1132±248 


5451±177 


Average RMSF of C" [nm] 


0.19±0.05 


0.14±0.06 


Radius of gyration of C° [nm] 


0.78±0.07 


0.83±0.03 


Lifetime of petpide-water H-bond [ps] 


0.87±0.05 


0.53±0.02 



TABLE I. Molecular Dynamics (MD) simulation results for 
an elastin mimetic peptide (VPGVG)3 under relaxed and 
strained states, as described in the text. Peptide-water El- 
bonds were counted when the donor-acceptor cutoff distance 
of 0.35nm and the hydrogen-donor-acceptor angle of 30° was 
satisfied. 



4ns MD simulations. The observed increase in backbone 
motion is also consistent with resulting MD simulations 
of strained (VPGVG)i8 [S]- In addition, the observed 
reduction of the amplitude of the backbone motion re- 



sulting from the strain, as quantified by the Root Mean 
Square Fluctuation (RMSF) of the C" shown in Table 1, 
is also consistent with the expected change according to 
the LEM model. 

Referring to Fig. 2a, the measured decrease in corre- 
lation time indicates an increase in motion and hence an 
increase in entropy of the proximal waters of hydration 
with increasing strain. This observation appears to be 
correlated with the observation that the lifetime of hydro- 
gen bonds formed between the model peptide (VPGVG)3 
and water decreases upon extension, as noted in Table 
1. Given that the deformation is applied to the protein 
backbone, it is clear that the changes to the dynamics of 
the proximal waters of hydration are driven by the in- 
crease in the frequency of libration of the backbone. The 
increase in the entropy of the proximal waters of hydra- 
tion and reduction in entropy of the backbone upon me- 
chanical strain therefore suggests that the driving mech- 
anism of elasticity is a change in the backbone entropy. 
While the experimental data suggest that the proximal 
waters of hydration in elastin increase mobility upon in- 
creasing strain, the presence of water is a required, albeit 
complex element for the elasticity of elastin. Ultimately 
the solvent gives rise to backbone mobility whose increase 
in librational frequency with increasing strain reduces the 
backbone entropy and drives the system back to a relaxed 
state. 
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